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Abstract. The gravitational wave emission by a distorted 
rotating fluid star is computed. The distortion is supposed 
to be symmetric around some axis inclined with respect 
to the rotation axis. In the general case, the gravitational 
radiation is emitted at two frequencies: Q and 257, where 
is the rotation frequency. The obtained formulae are ap- 
plied to the specific case of a neutron star distorted by 
its own magnetic field. Assuming that the period deriva- 
tive P of pulsars is a measure of their magnetic dipole 
moment, the gravitational wave amplitude can be related 
to the observable parameters P and P and to a factor /3 
which measures the efficiency of a given magnetic dipole 
moment in distorting the star. /3 depends on the nuclear 
matter equation of state and on the magnetic field distri- 
bution. The amplitude at the frequency 2^1, expressed in 
terms of P, P and /3, is independent of the angle a be- 
tween the magnetic axis and the rotation axis, whereas at 
the frequency Q, the amplitude increases as a decreases. 
The value of f3 for specific models of magnetic field distri- 
butions has been computed by means of a numerical code 
giving self-consistent models of magnetized neutron stars 
within general relativity. It is found that the distortion 
at fixed magnetic dipole moment is very dependent of the 
magnetic field distribution; a stochastic magnetic field or a 
superconductor stellar interior greatly increases (3 with re- 
spect to the uniformly magnetized perfect conductor case 
and might lead to gravitational waves detectable by the 
VIRGO or LIGO interferometers. The amplitude modu- 
lation of the signal induced by the daily rotation of the 
Earth has been computed and specified to the case of the 
Crab pulsar and VIRGO detector. 

Key words: gravitation - magnetic fields - gravitational 
radiation - stars: neutron - pulsars: general - numerical 
methods: spectral 
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1. Introduction 

Rapidly rotating neutron stars (pulsars) might be an im- 
portant source of continuous gravitational waves in the fre- 
quency bandwidth of the forthcoming LIGO and VIRGO 
interferometric detectors (cf. Bonazzola & Marck 1994 for 
a recent review about the astrophysical sources these de- 
tectors may observe). It is well known that a stationary 
rotating body, perfectly symmetric with respect to its ro- 
tation axis does not emit any gravitational wave. Thus 
in order to radiate gravitationally a pulsar must deviate 
from axisymmetry. Various kinds of pulsar asymmetries 
have been suggested in the literature: first the crust of a 
neutron star is solid, so that its shape may not be nec- 
essarily axisymmetric under the effect of rotation, as it 
would be for a fluid: deviations from axisymmetry are 
supported by anisotropic stresses in the solid. The shape 
of the crust depends not only on the geological history 
of the neutron star, especially on the episode of crystal- 
lization of the crust, but also on star quakes. Due to its 
violent formation (supernova) or due to its environment 
(accretion disk), the rotation axis may not coincide with 
a principal axis of the neutron star moment of inertia and 
the star may precess (Pines & Shaham 1974). Even if it 
keeps a perfectly axisymmetric shape, a freely precessing 
body radiates gravitational waves (Zimmermann & Sze- 
denits 1979). Neutron stars are known to have important 
magnetic fields; magnetic pressure (Lorentz forces exerted 
on the conducting matter) can distort the star if the mag- 
netic axis is not aligned with the rotation axis, which is 
widely supposed to occur in order to explain the pulsar 
phenomenon. Another mechanism for producing asymme- 
tries is the development of non-axisymmetric instabilities 
in rapidly rotating neutron stars driven by the gravita- 
tional radiation reaction (CFS instablity, cf. e.g. Schutz 
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1987) or by nuclear matter viscosity (cf. e.g. Bonazzola, 
Frieben & Gourgoulhon 1996 and references therein). 

In the seventies it was widely thought that the solid 
part of neutron stars was large (see Sect. Ill of Goldreich 
1970). A rotating rigid body whose rotation axis does not 
coincide with a principal axis of the moment of inertia 
must have a precessional motion. Therefore, there have 
been a number of studies about the precession of neutron 
stars and the resulting gravitational wave emission (Zim- 
mermann 1978a, Zimmermann & Szedenits 1979, Zimmer- 
mann 1980, Alpar & Pines 1985, Barone et al. 1988, de 
Araujo et al. 1994). However, modern dense matter cal- 
culations reveal that neutron star interiors are completely 
liquid (see e.g. Haensel 1995). Since the crust represents 
only 1 to 2% of the stellar mass (Lorenz et al. 1993), it 
appears that most of the neutron star is in a liquid phase. 
In this case, the star is not expected to precess, or at least 
its precession frequency must be reduced by a factor 10^ 
as compared with the solid case (Pines & Shaham 1974). 

The gravitational radiation from a rotating fluid star 
has not been studied as much as that from a solid star. 
The only reference we are aware of is the work by Gal'tsov, 
Tsvetkov & Tsirulev (1984) and by Gal'tsov & Tsvetkov 
(1984) about the gravitational emission of a rotating New- 
tonian homegeneous fluid drop with an oblique magnetic 
field. These authors have shown that the gravitational 
waves are emitted at two frequencies: the rotation fre- 
quency and twice the rotation frequency. They did not give 
explicit formulae for the two gravitational wave amplitude 
h+ and hx but have instead computed the response of 
a heterodyne detector (consisting in a rotating dumbbell 
orientated at right angles to the direction of propagation 
of the wave) directly from the Riemann tensor. 

In the present article, we compute the gravitational ra- 
diation from a rotating fluid star, distorted along a certain 
direction by some process (for example an internal mag- 
netic field), the angle between this direction and the rota- 
tion axis being arbitrary (Sect. ||). We do not assume that 
the star is a Newtonian object: it can have a large gravita- 
tional field described by the theory of general relativity — 
which is relevant for neutron stars. As an illustration we 
give the specific example of a Newtonian incompressible 
fluid with a uniform internal mag netic field (Sect. |). We 
consider also the (more general) case where the deforma- 
tion is due to an internal magnetic field which is also re- 
sponsible for the P of the pulsar (electromagnetic braking) 
(Sect. §). Using the numerical code presented in (Bocquet, 
Bonazzola, Gourgoulhon & Novak 1995) for magnetized 
rotating neutron stars in general relativity, we compute 
the gravitational emission for a specific 1.4 Mq neutron 
star model and for various configurations of the internal 
magnetic field. Finally, we derive the response of an inter- 
ferometric detector to the gravitational signal, taking into 



account the daily change induced by the Earth rotation 
in the relative orientation of the detector arms and the 
source (Sect. ^). 

2. Generation of gravitational waves by a rotating 
fluid star 

2.1. Thome's quadrupole moment 

Neutron stars are fully relativistic objects, so that the 
standard quadrupole formula for gravitational wave gen- 
eration — which is derived for non relativistic sources (cf. 
Sect. 36.10 of Misner, Thorne & Wheeler 1973, hereafter 
MTW) — does not a priori apply. However Ipser (1971) 
has shown that the gravitational radiation from a slowly 
rotating and fully relativistic star is given by a formula 
structurally identical to the weak field standard formula, 
provided that the involved multipole moments are defined 
in an appropriate way. 

The leading term in th£ gravitational radiation field is 
then given by the formulaH: 
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where r is the distance to the source and Pij := Sij — 
riVj /r^ is the tensor of projection transverse to the line of 
sight. Eq. (|lj) holds for highly relativistic sources provided 
that Tij is the mass quadrupole moment defined as the 
quadrupolar part of the 1 /r^ term of the l/r expansion of 
the metric coefficient goo in an asymptotically Cartesian 
and mass centered (ACMC) coordinate system (Thorne 
1980). This latter is a very broad coordinate class, which 
includes the harmonic coordinates. The precise definition 
of Tij is the following one. The space around a neutron 
star of mass M, "typical" radius R and angular velocity 
f2 can be divided in three regions: 

— the strong-field region: r < a few GMjt?; 

— the weak-field near zone: max(i?, a few GM /c^) ^ r < 

c/n 

— the wave zone : r > c/fl 

In the wave zone, retardation effects are important 
whereas in the weak-field near zone, the gravitational 
field can be considered as quasi-stationary. Note that 
for the fastest millisecond pulsar to date, PSR 1937-f21, 
c/il ^ 80 km, whereas i? ~ 10 km and if its mass is 
M = 1.4Mq, GM/c^ ^ 2 km. Hence, even in this extreme 
case, the weak-field near zone is well defined. A coordinate 
system {t ) is said to be asymptotically Cartesian 

and mass centered (ACMC) to order 1 if the metric admits 
the following l/r expansion (r :— XiX^) in the weak-field 



^ the results presented in this reference are erroneous and are 
corrected in Zimmermann & Szedenits (1979) 



^ The following conventions are used; Latin indices 
k,l, . . .) range from 1 to 3 and a summation is to be per- 
formed on repeated indices. 
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near zone (Sect. XI of Thornc 1980): 
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where ai, a2, cesij, oiiij, and f32ijk are some constants. 
The coefficients M and J*^ in the above expansion are the 
star's total mass and angular momentum. They are inde- 
pendent of the (a;') coordinate choice, provided that they 
reduce to usual Cartesian coordinates at the asymptoti- 
cally flat inflnity. On the contrary, the coefficient Tij is not 
invariant under change of asymptotically Cartesian coor- 
dinates: it is invariant only with respect to the sub-class 
of ACMC coordinates. This quantity is called the mass 
quadrupole moment of the star. As discussed in Sect. XI 
of Thorne (1980) (cf. also Thorne & Giirsel 1983), is 
a flat-space-type tensor which "resides" in the weak-fleld 
near zone, which means that it can be manipulated as a 
tensor in flat space. This tensor is symmetric and trace- 
free. If the star has a gravitational field weak enough to 
be correctly described by the Newtonian theory of grav- 
itation, Tij is expressible as minus the trace-free part of 
the moment of inertia tensor /,,-: 
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lij being given by 
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In particular, coincides at the Newtonian limit with the 
quantity introduced in MTW. For highly relativistic 
configurations, can no longer be expressed as some 
integral over the star, as in Eqs. (||)-(^; it must be read 
from the 1 /r expansion of the metric coefficient goo in an 
ACMC coordinate system. In Appendix we give the 
transformation from a quasi-isotropic coordinate system, 
usually used in relativistic studies of rotating neutron stars 
(cf the discussion in Sect 2 of Bonazzola, Gourgoulhon, 
Salgado & Marck 1993), to an ACMC coordinate system. 

2.2. Slightly deformed rotating star 

Very tight limits have been set on the deformation of 
neutron stars by pulsar timings: the relative deviation 
from axisymmetry is at most 10"'^ (cf. Table 1 of New 
et al. 1995). This number is certainly overestimated by 
several orders of magnitude for it is obtained by assum- 
ing that all the observed period derivative P is due to 



the gravitational emission, whereas it is much more likely 
to be due to the electromagnetic emission. Indeed, the 
present measurements of the braking index n of pulsars are 
n = 2.509 ± 0.001 for the Crab, n = 2.01 ± 0.02 for PSR 
0540-69 and n = 2.837 ± 0.001 for PSR 1509-58 (cf. Mus- 
limov & Page 1996 for references). Thus they are closer 
to n = 3 (magnetic dipole radiation) than io n — b (grav- 
itational radiation, cf. Table 9.1 of Manchester & Taylor 
1977). 

For such small deformations, we make the assumption 
that the quadrupole moment Tij can be linearly decom- 
posed into the sum of two pieces: 
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Tlj^ is the quadrupole moment due to rotation: X[°* = if 
the configuration is static. Tfj^^ is the quadrupole moment 
due to some process that distorts the star, for example 
an internal magnetic field or anisotropic stresses for the 
nuclear interactions. 

The fact that the star does not precess implies 
that there exists some ACMC coordinate system = 
{t,x,y,z) such that the components in this coordi- 
nate system do not depend upon x° = t and take the 
diagonal form 
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The line a; = y = is the rotation axis. For a small angular 
velocity fJ, X™* is a quadratic function of il. 

Let us make the following assumptions about the dis- 
tortion process: 

1. there exists some privileged direction, i.e. two of the 
three eigenvalues of Tfp^ are equal, i.e. there exists 
in the weak-field near zone a coordinate system (not 
necessarily ACMC) — {t, x, y, z) such that the com- 
ponents T^-^^ are 
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2. the deformation is rotating with the star, i.e. the weak- 
field near zone vector associated with the coordinate 
z is rotating at the angular velocity 17 at a fixed angle 
OL from the weak-field near zone vector associated 
with the coordinate z (e^ is along the rotation axis) 
(cf. Fig. I): 

Gi = cosaGz + sina [cos(p(f) ea; -I- sin</?(t) Gj,] , (10) 
with 

^(i) = f}(t - to) . (11) 
In particular, the above assumptions are satisfied by a 
magnetic field symmetric with respect to some axis: g^ is 
then parallel to the magnetic dipole moment vector A4. 
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taken. One obtains, using Eq. (O 



„ /2sinacos2iy9(t) 
-Xlf'rj^sinck 2sinasin2^(t) 
\ cosQ!sin95(i) 

2sinasin293(t) cosasin(y9(t) 
— 2sinacos2(/5(f) — cos a cos (p{t) 
— cos a cos ip{t) 



(15) 



Fig. 1. Geometry of the distorted neutron star. 



The next step consists in taking the transverse traceless 
projection of I^'^"'*. Let us denote by i the angle between 
the neutron star's rotation axis and the direction n 
from the star's centre to the Earth {i is called the line of 
sight inclination) (cf. Fig. ||). Without any loss of general- 
ity, we can choose the ACMC coordinates (t, a;, y, z) such 
that n lies in the (ej,,ez) plane. The components with re- 
spect to {t,x,y,z) of the transverse projection operator 
Pa = (5,o — n,no are then 



2. 3. Application of the quadrupole formula 



ing the time constancy of X[™ , we get 



Inserting Eq. into the quadrupole formula (|^) and us 

('-;) 



XT 2(j 1 
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T-dist 



(12) 



In order to apply this formula, one must first express 
Tfj^^ in the ACMC coordinates (t, x, y, z), from its compo- 
nents (|^) in the coordinates (i, x,y, z). The transformation 
{t,x,y,z) — > (t,x,y,z) is simply the composition of two 
rotations (cf . Fig. |l|) : one of angle a around an axis e^' 
in the {ex:,ey) plane and rotating with the star, and one 
of angle ip{t) = Q{t — tQ) around e^. The transformation 
matrix is then: 



^cos ip{t) ~ sirup (t) 
P= I sin ip{t) cos ip{t) 
1 





cos a sin a 
— sin a cos a 



.(13) 



The tensor transformation law writes X' 
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*P, where I"^"^* is the matrix formed by the components 
Zj^* and jdist jg |-]-^g rnatrix formed by the components 

■^^dist 



P — 



cos^ i 



> — sin i cos i 
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From Eqs. ( p^ and (15), the computation of the right- 
hand side of the quadrupole formula ([l^) is straightfor- 
ward, though somewhat tedious. The result is 



with 



h+e+ + h^ef^ 



'10 
— cos^ i sin i cos i 
J sin i cos i 

cos i — sm I 
cos i 
— sin I 
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(18) 
(19) 



and 

h. = 



ho sin a 



■ cos a sini cosi cos fl{t — to) 



sma- 



1 + cos i 



cos 2il{t — to) 



[Eq. 
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Performing this matrix product leads to = ho sin a - cos a sin i sin SI (i - to) 



z; 



■dist 
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3sin^ asin^ ip{t) - 1 



sin a cosi sin 2Q{t — to] 



= -X? 
2 



3/2 sin^ asin2(p{t) 
3 sin a cos a sin ip{t) 

—3/2 sin^ asin2ip{t) — 3sinacosasin</5(t) \ 
Ssin^ acos^ (p(t) — 1 3 sin a cos a cos (/?(<) 
3 sina cos a cos iy9(t) Scos^a — 1 J 



(14) 



where 

ho := - 



6G distill 
c4 r 



(20) 



(21) 



(22) 



In order to apply the quadrupole formula, the second 
derivative with respect to t of this expression must be 



Note that in the above expressions, Eq. (^l]) has been sub- 
stituted for ip[t) and that the retardation term — r/c has 
been incorporated into the constant to. 
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8.4- Discussion 

From the formulae (20)-(21), it is clear that there is no 
gravitational emission if the distortion axis is aligned with 
the rotation axis (a = or tt) . If both axes are perpendic- 
ular {a = 7r/2), the gravitational emission is monochro- 
matic at twice the rotation frequency. In the general case 
(0 < \a\ < 7r/2), it contains two frequencies: and 251. 
For small values of a the emission at f2 is dominant. 

It may be noticed that Eqs. (20)-(p^ are structurally 
equivalent to Eq. (1) of Zimmermann & Szedenits (1979) 
(hereafter ZS), although this latter work is based on a 
different physical hypothesis (Newtonian precessing rigid 
star). In order to compare precisely Eqs. (20)-(p2|) and 
Eq. (1) of ZS, some slight re-arrangements must be per- 
formed. First, the ellipticity e defined by ZS is linked to 
our Izz by he = —3/2122 [cf. Eq. (||)], this explains the 
factor 2 in front of the right-hand side of Eq. (1) of ZS in- 
stead of the factor —6 in Eq. (^). Other apparent differ- 
ences are actually due to different conventions: the origin 
of time in ZS is the instant when the equivalent of the de- 
formation axis is at its farthest position from the observer, 
whereas in our case it corresponds to its nearest position. 
One must then make t —^ t+i: in order to compare the two 

formulae. Moreover the choice for the matrices eti and ef,- 

''J ''J 

are exactly opposite in both approaches [cf. Sect. III.B of 
Zimmermann (1980)], so that the transforms /i+ — > — /i+ 
and /ix ^ must also be performed. When all this 

is done, Eqs. (20)-(|2|) and Eq. (1) of ZS appear to have 
exactly the same structure. However, the physical signifi- 
cance is different: the frequency u> which appears in Eq. (1) 
of ZS differs from the pulsar frequency by the body-frame 
precessional frequency whereas in Eqs. (20)-(|2^), Vl is ex- 
actly the pulsar frequency. The angle 9 of ZS, which in 
their Eq. (1) takes the place of our angle a in Eqs. (20)- 
(p2|), is the angle between the total angular momentum J 
and the star's third principal axis, whereas our a is the 
angle between the rotation axis and the direction of the 
distortion, which even in the Newtonian case, does not co- 
incide with any of the principal axis of the body (except 
in the non-rotating case). 

As special cases of Eqs. (20)-(p2|), one may recover re- 
sults previously published in the literature. For instance, 
the case of a triaxial star rotating about a principal axis 
of its moment of inertia tensor can be obtained by set- 
ting a = 7r/2 in Eqs. (20)-(p2|). The result can be com- 
pared with Eqs. (48) and (54) of Thorne (1987), noticing 
that / in these equations is fi/Tr and that X^x and Xyy 
of Thorne (1987) are related to our In by I^^ = — Iii/2 
and lyy = Izi ■ The two formulae appear then to be identi- 
cal, as expected. If in addition to a = 7r/2, the inclination 
angle i is set to zero, Eqs. (20)-(|2^) reduce to the for- 
mula used in the recent work by New et al. (1995) [their 
Eq. (5)]. Note that in the two studies mentionned above, 
the gravitational waves are emitted at the frequency 2U, 
only, due to a = 7r/2 or equivalently due to the fact that 



the rotation axis coincides with a principal axis of the 
moment of inertia tensor. Let us stress again that in our 
(more general) case, the gravitational radiation contains 
two frequencies: J7 and 2Vl. 

2.5. Numerical estimates 

In order to describe the star deformation by a dimension- 
less quantity, let us introduce instead of I'lf'^ ellipticity 
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where / is the moment of inertia with respect to the ro- 
tation axis, defined as 



/ := J /VI 



(24) 



J being the star angular momentum. The definition (24) 
is valid even in highly relativistic cases, provided that the 
star distortion is small — as we suppose throughout this 
work. Indeed, in this case the configuration is essentially 
axisymmetric and the angular momentum J is well de- 
fined (cf. e.g. the discussion in Wald (1984), p. 297). The 
factor —3/2 in Eq. ( p3| ) is introduced in order to recover 
the classical definition of the ellipticity at the Newtonian 
limit (see e.g. Shapiro & Teukolsky 1983). 

Let us introduce also the rotation period P = 2w/i^, 
since observational data about pulsars are usually pre- 
sented with P instead of 17. 

Inserting Eq. ( p3| ) in expression (|2^) for the character- 
istic gravitational wave amplitude leads to 



ho 



P^r 



(25) 



Replacing the physical constants by their numerical values 
results in 
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Note that I = 10^* kgm^ is a representative value for the 
moment of inertia of a 1.4 neutron star [see Fig. 12 of 
Arnett & Bowers (1977)]. 

For the Crab pulsar, P = 33 ms and r = 2 kpc, so that 
Eq. (Eq) becomes 



(27) 



/i^''^'' = 1.89 X 10^27 



ID- 



IOMS kgm^ 

For the Vela pulsar, P = 89 ms and r = 0.5 kpc, hence 



I, Vela 
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(28) 
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For the millisecond pulsarB PSR 1957+20, P = 1.61 ms 
and r ~ 1.5 kpc, hence 



^1057+20 ^ Qg ^ -^Q- 



24 


r J 1 




e 








[lO-eJ 



(29) 



At first glance, PSR 1957+20 seems to be a much more 
favorable candidate than the Crab or Vela. However, in 
the above formula, e is in units of 10~^ and the very low 
value of the period derivative P of PSR 1957+20 implies 
that its e is at most 1.6 x 10~^ (New et al. 1995). Hence 
the maximum amplitude one can expect for this pulsar is 
^1957+20 _ X ^Q-n j^o^ 08 X 10-24 as Eq. (|9|) 

might suggest. 



3. The incompressible magnetized fluid example 

As an illustration, let us consider the specific example 
taken by Gal'tsov, Tsvetkov & Tsirulev (1984) (hereafter 
GTT) and Gal'tsov & Tsevtkov (1984). In their work, a 
pulsar is idealized as a rigidly rotating Newtonian body 
made of an incompressible fluid and endoved with a mag- 
netic field which is uniform inside the star and dipolar 
oustide it, the magnetic dipole moment being inclined by 
an angle a with respect to the rotation axis. 

The rotation rate is supposed to be far from the mass- 
shedding limit so that the departure from spherical sym- 
metry is small. Moreover, the magnetic energy is assumed 
to be much lower than the rotational kinetic energy, which 
is satisfied by realistic configurations. Under these hy- 
potheses, the star takes the shape of a (quasi-spherical) 
triaxial ellipsoid. The gravitational potential can be then 
calculated from the analytical formulae of Chandrasekhar 
(1969). The shape of the ellipsoid is deduced from the first 
integral of the equation of motion (including the magnetic 
pressure) and the magnetic field matching conditions at 
the stellar surface. It is found that the ellipsoid is deter- 
mined by the equation + aij)X'^X^ = R^, where (i) 
R is the mean radius of the (quasi-spherical) ellipsoid, (ii) 

are Cartesian coordinates in a co-moving frame, i.e. ro- 
tating at the angular velocity witf^ respect to an inertial 
frame and (iii) the aij are given hy 



15 n.n^ 



45^0 MiMj 



327r2 i?8pwj 



(30) 



In this expression, p is the constant mass density of the 
star, wj = ^/AnGp is the Jeans frequency, Qi and Mi are 
respectively the components of the angular velocity vec- 
tor and the components of the magnetic dipole moment 

We do not consider the "historical" millisecond pulsar PSR 
1937+21 for it is more than twice farther away. 
* Note that the value of atj presented in GTT [their 
Eq. (2.15)] is erroneous: the fli part has the wrong sign and 
the Mi part should contain a factor instead of the i?*. 
This latter error has been corrected in Gal'tsov & Tsvetkov 
(1984), but not the former one. 



with respect to the X' coordinates: fii — (0, 0, H.) and 
Mi ~ {0, M sina, M cosa). By diagonalizing the matrix 
aij given by Eq. ( |30| ) , one obtains the principal axes of the 
ellipsoid and the values of the three semi-axes, ai, a2 and 
03. From these quantities the moment of inertia tensor 
lij can be computed by evaluating the integral (H) in the 
frame of the principal axes. The quadrupole moment 1^ is 
then immediately deduced via Eq. (^). Transforming the 
result in the inertial frame leads to the form (^ of lij , in- 
cidently demonstrating this formula in the particular case 
under consideration, with the form (^) for X™* with 



jrot _ " 

~ 3G 

and with the form (14) for I^^^^ with 
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dist 



poM'^ 



16n^GpR^ 



(31) 



(32) 



The moment of inertia with respect to the rotation axis of 
the homogeneous star is / = SnpR^ so that Eq. (23) 
gives the ellipticity: 



45 



R2 

^polo 



647r fioGp^R^ 



(33) 



Note that in this formula, the magnetic dipole moment M 
has been expressed in terms of the North pole magnetic 
field Bpoic = /io/(47r) 27W/i?3. 

For numerical estimates, let us take typical values for 
neutron stars: Spoio = 10^ T, M = 1.4 Mq and i? = 10 
km. Then e ~ 6.0 x 10^^^. This is a very tiny value, 
which leads to a gravitational wave amplitude of only 
/lo ~ 3 X 10-30 for P = 10 ms and r = 1 kpc [cf. Eq. (|6|)]. 
However the above model is a very simplified one. It can 
be expected that relaxing the assumptions of (i) incom- 
pressible fluid, (ii) Newtonian gravity and (iii) uniform 
internal magnetic fleld, may lead to a greater value of e. 

4. Magnetic field induced deformation 

Emission formula 

The situation considered in the preceding section is a very 
simplified one. However, one may consider that the ob- 
tained form of the deformation, Eq. (|3^), is qualitatively 
the same for realistic magnetized neutron star models, 
i.e. a compressible perfect fluid obeying a "sophisticated" 
equation of state resulting from nuclear physics calcula- 
tions, and involving general relativity. More precisely, we 
consider that the magnetic fleld induced deformation is a 
quadratic function of the amplitude of the magnetic dipole 
moment, M, as in Eq. (B2h: 



,M^ 
Mf 



h 42 ^ 



(34) 
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where Mo has the dimension of a magnetic dipole mo- 
ment in order to make the coefficient /? dimcnsionless. Let 
us choose Mq = (An/fio) GP/R"^ where R is the circum- 
ferential equatorial radius of the starcl. Hence 



Ho M^R^ 
An GP 



(35) 



Provided the magnetic field amplitude does not take (un- 
realistic) huge values (> 10^'* T), this formula is certainly 
true, even if the magnetic field structure is quite compli- 
cated, depending upon the assumed electromagnetic prop- 
erties of the fiuid: normal conductor, superconductor, fer- 
romagnetic... The coefficient /? measures the efficiency of 
this magnetic structure in distorting the star. In the fol- 
lowing, we shall call /? the magnetic distortion factor. For 
the simplified model considered in Sect. || (incompressible 
fluid, uniform internal magnetic field), (3 = 1/5. 

As argued in Sect. 2.2, the observed spin down of radio 
pulsars is very certainly due to the low frequency magnetic 
dipole radiation. M is then linked to the observed pulsar 
period P and period derivative P by [cf e.g. Eq. (6.10.26) 
of Straumann (1984)] 



-Bpoic = 10^ T and P = 10 ms has a moment of inertia 
/ = 1.1 X 10'^* kgni^, a factor f3 — 1/5, and magnetic 
dipole moment = 5.0 x 10^'^ Am^. From Eq. (56|), the 
corresponding period derivative is P = 2.2 x lO^^^in^ a. 
Equation ( |4C| ) gives then /iq = 3 x 10"'^° for r = 1 kpc, in 
agreement with the result obtained in Sect. ^. 

The main feature of the emission formulae (38)- (39) is 
that the amplitude of the gravitational radiation at the 
frequency 2Q does not depend upon the (unknown) angle 
a. 

4-.2. Discussion 

Among the 706 pulsars of the catalog by Taylor et al. 
(1995, 1993), the highest value of Hq at fixed a, (3 and 
R, as given by Eq. ([l0|), is achieved by the Crab pulsar 
(P = 33 ms, P = 4.21 x lO'^^^ r = 2 kpc), followed 
by Vela (P = 89 ms, P = 1.25 x lO^^^ r = 0.5 kpc) 
and PSR 1509-58 (P ^ 151 ms, P = 1.54 x lO^^^^ r = 
4.4 kpc): 



An 3c-^ IPP 



^ sin^ a 



(36) 



where a is the angle between the magnetic dipole moment 
and the rotation axis. For highly relativistic configura- 
tions, the vector A4 is defined in the weak-field near zone 
(cf. Sect. 2.5 of Bocquet et al. 1995), so is a. Inserting 
Eqs 



I^Crab 
"■0 



,Vcla 



I 1509-58 

"■0 



t 1957+20 

"o 



4.08 X 10--^^ 
1.81 X IQ-^^ 

1.50 X IQ-^^ 

4.51 X 10-3"^ 



R 



sin^ a 



10 km 

R P (3 



SVC? a 



10 km 

R P P 



sin^ a 



10 km 



10 km. 



sin^ a 



(41) 
(42) 
(43) 
(44) 



35|) and (|36|) into Eq, 
wave amplitude 



5[) leads to the gravitational We have added to the list the millisecond pulsar PSR 



ho = 613- 



R^P 



crP sin a 
so that formulas (20) and (21) become 

P2p 



(37) 



h+ = 6/3 



6/3 



crP 



sm I cos I 
2 tan a 

1 + cos^ i 



cos fl(t — to) 

cos 2n,{t — to) 



R^P\ sini 



crP 



sin il{t — to) 



2 tan a 
~ cos z sin 2fl{t — to) 



(38) 



(39) 



1957-1-20 (P = 1.61 ms, P = 1.68 x 10-^°, r = 1.5 kpc) 
considered in Sect. 2.5. From the above values, it appears 
that PSR 1957-1-20 is not a good candidate. This is not 
suprising since it has a small magnetic field (yielding a 
low P). Even for the Crab and Vela pulsars, which have 
a large P, the ho values as given by Eqs. (41), (42) are, 
at first glance, not very encouraging. Let us recall that 
with the 10-22 Ez'^P expected sensitivity of the VIRGO 
experiment at the 30 Hz frequency (Bondu 1996; see also 
Fig. 9 of Bonazzola & Marck 1994), the minimal ampli- 
tude detectable within three years of integration is 



10 



-26 



(45) 



Equation (|3) 
form: 

ho = 6.48x10" 



can be cast in a numerically convenient 



.(40) 



30 P 


\ ^ 1 


2 


"kpc 




ms' 




\ ^ 1 


■ 2 

sm a 


aOkmJ 




. r . 




I P \ 




ilO-isJ 



As a check of this equation, let us consider again the sim- 
plified case of Sect. The uniformly magnetized homo- 
geneous Newtonian star with M = 1.4 M©, P = 10 km. 



Comparing this number with Eqs. (38)-(39) and (41)-(42), 
one realizes that in order to lead to a detectable signal, the 
angle a must be small and/or the distortion factor /3 must 
be large. In the former case, the emission is mainly at the 
frequency fi. From Eqs. (38)-(39), the gravitational wave 
amplitude can even be arbitrary large if a ^ 0. However, 
if a is too small, let say a < 10-^, the simple magnetic 
braking formula (36) certainly breaks down. So oiie can- 
not rely on a tiny a to yield a detectable amplitudea. The 



® For Newtonian stars, R is simply the equatorial radius; for 
relativistic stars, R is the length of the equator, as measured 
by a locally non rotating observer, divided by 2tt. 



® From the observed pulse profile and polarization of pulsars, 
values of a can be inferred; they spread all the range between 
and 7r/2 (Lyne & Manchester 1988, Rankin 1990) 
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alternative solution is to have a large p. Let us recall that 
for an incompressible fluid with a uniform magnetic field, 
(3 — 1/5 (Sect. H). In the following section, we give the 
(3 coefficients computed for more realistic models (com- 
pressible fluid, realistic equation of state, general relativ- 
ity taken into account) with various distribution of the 
magnetic field. 



Magnetic field lines 




X [ km ] 

Fig. 2. Magnetic field lines generated by the current distri- 
bution corresponding to the choice / = const. The thick line 
denotes the star's surface. The distortion factor of this config- 
uration is /3 = 1.01. 



4-. 3. Numerical results 

We have developed a numerical code to compute the de- 
formation of magnetized neutron stars within general rel- 
ativity. This code is an extension of that presented in Boc- 
quet et al. (1995) (hereafter BBGN). The main improve- 
ments are (i) the use of an arbitrary number of grids to de- 
scribe the stellar interior, which allows a greater diversity 
of magnetic field configurations, and (ii) the possibility of 
a type I superconductor interior. We report to BBGN for 
details about the relativistic formulation of Maxwell equa- 
tions and the technique to solve them. Let us simply recall 
here that the obtained solutions are fully relativistic and 
self-consistent, all the effects of the electromagnetic field 
on the star's equilibrium (Lorentz force, spacetime curva- 
ture generated by the electromagnetic stress-energy) be- 
ing taken into account. The magnetic field is axisymmetric 
and poloidal. The numerical technique is based on a spec- 
tral method (numerical details can be found in Bonazzola 
et al. 1993). 



Magnetic field lines 




X [ km ] 

Fig. 3. Magnetic field lines generated by a current distribu- 
tion localized in the core of the star and corresponding to the 
choice / = const. The thick line denotes the star's surface and 
the dashed line the external limit of the electric current distri- 
bution. The distortion factor corresponding to this situation is 
P = 0.70. 

Thanks to the splitting (^, we do not need to take 
into account the rotation to compute the Xfj^*^ induced 
by the magnetic field. Consequently we consider static 
magnetized neutron star models. The reference (non- 
magnetized) configuration is taken to be a IAMq static 
neutron star built with the equation of state UV14 + TNI 
of Wiringa, Fiks & Fabricini (1988). This latter is a mod- 
ern and medium stiffness equation of state (cf. Sect. 4.1.2 
of Salgado et al. 1994). The circumferential radius is 
R = 10.92 km, the baryon mass 1.56 Mq, the moment 
of inertia / — 1.23 x 10"^^ kgm^ and the central value of 
goo is 0.36, which shows that such an object is highly rel- 
ativistic. Various magnetic field configurations have been 
considered; the most representative of which are presented 
hereafter. 

4.3.1. Normal case 

Let us first consider the case of a perfectly conducting 
interior (normal matter, non-superconducting). As dis- 
cussed in BBGN, the electric current distribution j can- 
not be arbitrary in order to lead to a stationary con- 
figuration: it must be related to the covariant (p com- 
ponent of the electromagnetic potential vector, A^, by 
j"^ — rij* = (e + p)f{A^), where e and p are respectively 
the proper energy density and pressure of the fluid and 
/ is an arbitrary function. The simplest magnetic config- 
uration is given by the choice f{x) — const. It results in 
electric currents in the whole star with a maximum value 
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at half the steUar radius in the equatorial plane. The cor- 
responding magnetic field distribution is shown (in coor- 
dinate space) in Fig. ^. The resulting distortion factor is 
P = 1.01, which is above the 1/5 value of the uniform 
magnetic field/incompressible fluid Newtonian model con- 
sidered in Sect. ^ but still very low. This is not surprising 
since the magnetic field configuration has a very simple 
structure: it is certainly not the configuration which max- 
imizes the deformation at fixed magnetic dipole moment. 

Changing the function / does not lead to a dramatic 
increase in /3: for f{x) = ai(a;exp(— x^) -|- a2), we get 
P = 1.07 and for f{x) = ai cos{x/a2), P = 1.29. 

Our multi-grid code allows to study localized distri- 
bution of the electric current. Figure || corresponds to the 
electric current distribution given by / = const from r ~ 
up to — 0.5 Toq, where is the r coordinate of the 
equator, and to no electric current for r > r^,. The distor- 
tion factor is only /3 = 0.70. This can be understood since 
in the absence of electric current, there is no Lorentz force 
in the outer part of the star. For electric currents concen- 
trated deep in the stellar core, the situation is more favor- 
able. Indeed, since in the regions where j = the magnetic 
field B falls off as ~ r~^, a moderate value of can cor- 
respond to an important value of B at the stellar centre, 
leading to a substantial deformation of the core, that grav- 
itationally influences the rest of the star. For instance, for 

= 0.1 Teq, we get P = 5.86. 

The opposite situation corresponds to electric currents 
localized in the neutron star crust only. Figure ^ presents 
one such configuration: the electric current is limited to 
the zone r > r^, — 0.9 Toq {f{x) = const.). The resulting 
distortion factor is /3 = 8.84. 



4.3.2. Type I superconductor 

Let us consider the case of a superconducting interior, of 
type I, which means that all the magnetic field has been 
expulsed from the superconducting region. In the config- 
uration depicted in Fig. ||, the neutron star interior is su- 
perconducting up to = 0.9 Teq. For r > r,, the matter 
is assumed to be a perfect conductor carrying an electric 
current which corresponds to f{x) = const. The resulting 
distortion factor is much higher than in the normal case: 
P — 157. For = 0.95 roq, P is even higher: P = 517. 



4.3.3. Counter-rotating electric currents 

The above values of P, of the order 10^ — 10'^, though much 
higher than in the simple normal case (Sect. [4.3.1 ), are 
still too low to lead to an amplitude detectable by the first 
generation of interferometric detectors in the case of the 
Crab or Vela pulsar [cf. Eqs. (41), (42) and (^)]. It is clear 
that the more disordered the magnetic field the higher /?, 
the extreme situation being reached by a stochastic mag- 
netic field: the total magnetic dipole moment M almost 
vanishes, in agreement with the observed small value of P, 



Magnetic field lines 




Fig. 4. Magnetic field lines generated by a current distribu- 
tion localized in the crust of the star and corresponding to the 
choice / — const. The thick line denotes the star's surface and 
the dashed line the internal limit of the electric current distri- 
bution. The distortion factor corresponding to this situation is 
f3 = 8.84. 



Magnetic field lines 




Fig. 5. Magnetic field lines generated by a current distribution 
exterior to a type I superconducting core. The thick line de- 
notes the star's surface and the dashed line the external limit of 
the superconducting region. The distortion factor correspond- 
ing to this situation is /3 = 157. 
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Magnetic field lines 




-10 





km 



10 



Fig. 6. Magnetic field lines generated by some electric current 
rotating in one direction for < r < 0.3 req and in tiie opposite 
direction for 0.9 rcq < r < req. The thick line denotes the 
star's surface and the dashed lines the limits of the electric 
current distributions. The distortion factor corresponding to 
this situation is f) — 5.7 x 10^. 



whereas the mean value of throughout the star is huge. 
Note that, according to Thompson & Duncan (1993), tur- 
bulent dynamo amphfication driven by convection in the 
new-born neutron star may generate small scale magnetic 
fields as strong as 3 x 10^^ T with low values of -Bdipoic 
outside the star and hence a large (3. 

In order to mimic such a stochastic magnetic field, let 
us consider the case of counter-rotating electric currents, 
namely j is given by (i) f{x) = fi for < r < r^, where 
/i is a constant, (ii) f{x) — ji for r'^<r< rcq, where ji 
is a constant with the opposite sign than /i and (iii) no 
electric current for r* < r < rj . Figure ^ corresponds to 



the case r* = 0.3 To 



0.9 req and /1//2 = -5.45. The 



resulting distortion factor is /3 = 5.7 x 10'^. The value of 
magnetic dipole moment, M. = 5.1 x 10^^ Am^, is similar 
to that of the Crab pulsar (assuming a ~ 1), but the 
amplitude of the magnetic field at the star's centre, -Be ~ 
5.7 X 10^2 T = 5.7 X 10^6 G is - 10"* times higher than the 
polar value deduced from a simple dipole model. Clearly 
for such configurations /3 can be made arbitrary large by 
adjusting the parameters /i and /2. 

4.3.4. Type II superconductor 

It is not clear if the proto ns of the neutron star interior 
form a type I (Sect. 4.3.2| ) or a type II superconductor 
(P. Haensel, private communication). In the latter case, 
the magnetic field inside the star is organized in an array 
of quantized magnetic flux tubes, each tube containing a 



magnetic field ^ 10^^ T (Ruderman 1991). Besides, 
the neutrons constitute a superfluid, with quantized vor- 
tices. As the neutron star spins down, the neutron vortices 
migrate away from the rotation axis. As discussed by Rud- 
erman (1991, 1994), the magnetic flux tubes are forced to 
move with them. However, they are pinned in the highly 
conducting crust. This results in crustal stresses of the or- 
der Be5/2/io [Ruderman 1991, Eq. (10)], where B is the 
mean value of the magnetic field in the crust {B ~ 10® T 
for typical pulsars). This means that the crust is submit- 
ted to stresses ^ 10'^ higher than in the unifo rmly dis- 
tributed magnetic field considered in Sect. 4.3.1 (compare 
B^Bjl^iQ with B^/2^o)- The magnetic distortion factor /3 
should increase in the same proportion. We have not done 
any numerical computation to confirm this but plan to 
study type II superconducting interiors in a future paper. 

5. Signal received by an interferometric detector 

Due to the weakness of the expected gravitational signal 
from pulsars, long integration times, typically of the or- 
der of the year, are required to extract the signal out of 
the noise. For such a long observing time, the motion of 
the Earth must be taken into account for it modifies the 
position of the pulsar with respect to the antenna pat- 
tern of the interferometric detector. It results in a mod- 
ulation of both the amplitude and the frequency of the 
signal (Jotania, Valluri & Dhurandhar 1995). Note that if 
one is searching for a known pulsar, the frequency mod- 
ulation (Doppler shift) can be obtained readily from the 
radio observations. 

In this section, we examine the daily modulation of the 
signal amplitude due to the Earth's rotation. The special 
case of an interferometric detector situated on the equator 
with arms symmetrically placed about the North-South 
direction and a gravitational wave coming from the North- 
ern celestial pole has been treated by Jotania et al. (1995). 

5.1. Beam-pattern jactors 

Let (ej, e^, ej) be the orthonormal frame associated with 
the gravitational wave: e^ is perpendicular to the wave 
plane and parallel to the "line of sight" , being orientated 
from the neutron star centre to the Earth centre; ej is the 
unit vector of the wave plane which is perpendicular to the 
star's rotation axis. In the (e^, e^, ej) frame the transverse 
traceless gravitational wave is expressible from its two po- 
larization modes /14. and /ix as given by Eqs. (20)-(21) or 
(38)-(39): 



/i+ (ej ( 



- hy, (ej ( 



.(46) 



The response of an interferometric detector of the 
VIRGO/LIGO type to the above gravitational wave de- 
pends upon the relative orientation of wave frame with 
respect to the detector's arms. Let (i&Xt^y t^z) be the 
orthonormal frame such that e^ is perpendicular to the 



10 



S. Bonazzola & E. Gourgoulhon: Gravitational waves from pulsars 



11 




w 



Fig. 7. Relative orientation of various orthonormal frames introduced in the text: {ex,ey,ez) is the neutron star inertial 
(non-rotating) frame, associated with the ACMC coordinates {x,y,z); {e^jBy) define the plane S perpendicular to the rotation 
axis, {est, By, ej) is the frame associated with the gravitational wave received on Earth; (ej, ey) define the wave plane W. Note 
that 64 = 63;. (6x,6y,6z) is the interferometric detector frame; (6x,6y) define the detector plane D. 



detector plane, pointing toward the zenith and and 
By are unit vectors along the two detector's arms. Let 
(9, 'I') be the three Euler angles which specify the po- 
sition of the wave frame (ej,ey,ej) with respect to the 
detector frame {ex,eY,Gz) (cf. Fig. ^). The signal mea- 
sured by the detector is 

/i(i) = i^+(e,$,«');i+(i) + Fx(e,$,*)/ix(t) , (47) 

with the following beam-pattern factors [cf. e.g. 
Eqs. (103)-(104) of Thorne (1987) or Eq. (7a) of Dhurand- 
har & Tinto (1988)] : 

i^+(e,$,*) i(l + cos^e)cos2$cos2\E' 

-cose sin 2$ sin 2* (48) 

Fx(e,$,*) = i(l + cos2e)cos2<I>sin2* 

-I- cose sin 2$ cos 2 . (49) 

The Euler angles (e, $, ^) arc not constant because of the 
motion of the detector with respect to the source induced 
by the Earth's diurnal rotation and revolution around the 
Sun. To compute their variations let us introduce the "ce- 
lestial sphere frame" {ex> ,eY' ,ez') such that ez' is along 
the Earth's rotation axis, pointing toward the North pole, 
ex' and ey/ are in the Earth's equatorial plane, ex' point- 
ing toward the vernal point (i.e. ex' is along the inter- 
section of the Earth's equatorial plane with the Earth's 
orbital plane). For time scales of the order of the year, 
{ex' ,eY' ,ez') can be considered as fixed with respect to 



the (approximatively) inertial frame containing the solar 
system barycentre and the neutron star centre. The wave 
frame (65,6^,6^) is also fixed with respect to this iner- 
tial frame. Then let (6',$',^') be the Euler angles which 
specify the position of the wave frame with respect to the 
celestial sphere frame (cf. Fig. ||). e' and <&' are simply 
related to the equatorial coordinates of the Jieutron star 
on the celestial sphere (the right ascensiorfl a and the 
declination S) by (cf. Fig. ^): 

e' = ^ + I and $' = a - I . (50) 
The triad (e5,ej^,ej) is related to the triad (e^', ey;, e^/) 

by 

= A^j, ej, , (51) 

with the orthogonal matrix [cf. e.g. Eq. (4-46) of Goldstein 
(1980) with the substitution (|^)] 

(sin a cos '0 — cos a sin 5 sin ip 
— sin a sin -0 — cos a sin 5 cos 
— cos a cos 5 

— cos a cos "0 — sin a sin 5 sin ip cos 5 sin \ 
cos a sin -0 — sin a sin (5 cos '0 cos ^ cos ^ . (52) 
— sin a cos 5 — sin 6 I 

Let {bx" tGyii ,ez") be the orthonormal frame linked 
to the geographical location of the detector site, such that 

^ a bar is put on a to distinguish it from the angle a between 
the magnetic and rotation axis introduced earlier in the text. 
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Fig. 8. Relative orientation of the gravitational wave frame 
{es:,ey,ei) and the celestial sphere frame (e^^/, e^/, e^'). 



Fig. 9. Orientation of the detector arms with respect to the 
local North-South and West-East directions. The dashed line 
is the projection of the gravitational wave propagation direc- 
tion on the horizontal plane. The depicted configuration cor- 
responds to that of the VIRGO detector (A = -161° 33')- 



ez" = ez is the local vertical, ex" is in the North-South 
direction and ey is in the West-East direction (cf. Fig. 
The position of the wave frame [bx, ey,ez) with respect to 
the cardinal frame (ex" , ey, e^") is determined by the 
three Euler angles (O, where O and 5" are the same 

angles as those relative to the detector frame and entering 
Eqs. (48)-(49). The angle $" is related to $ by 



$ = + A , 



(53) 



where A is the azimuth of the X-arm of the detector, i.e. 
the angle between the South direction and the X-aim, as 
measured in the retrograd way. is directly related to 
the azimuth a of the source by 



TT 



(54) 



The minus sign which occurs in this relation comes from 
the fact that the azimuth is measured westwards from the 
South, hence in the inverse trigronometric way. Putting 
together relations (p3|) and (E^) we obtain 



$ = -a + A - 



(55) 



Let {6" , (p" , be the Euler angles which specify the po- 
sition of the cardinal frame {ex" ,eY" ,^z") with respect 
to the celestial sphere frame {ex' ,eY' ,^z')- One has im- 
mediatly tp" = — 7r/2 because ex" is orientated in the 
North-South direction. 6" is simply related to the latitude 
/ of the detector site by 



2 



(56) 



(p" is linked to the local sidereal time T (i.e. the angle 
between the local meridian and the vernal point) by 



(57) 



T can be expressed in terms of the sidereal time at Green- 
wich at Oh UT, rGrocn(O), and the local UT time t by 



T{t) =at + rGrco„(0) - L , 



(58) 



where cr = 1.00273790935 x 157h (Meeus 1991) is a con- 
version factor from mean solar time to sidereal time (hence 
(T accounts for the revolution of the Earth around the Sun) 
and L is the longitude of the detector site (the minus sign 
in front of L comes from the fact that the geographical 
longitudes are measured positively westwards). The triad 
{ex" jGy" ,ez") is related to the triad (ex', ey. e^') by 



with the orthogonal matrix 

(sin I cos T sin I sin T — cos I 
-sinT COST 
cos I COS T cos I sin T sin I 



(59) 



(60) 



Now the wave frame (ej,ey,ej) is related to the triad 
{ex",eY»,ez") by 



with the orthogonal matrix 

cos ^ cos $" — cos 9 sin sin \1/ 
C~^ j„ = I -sin^'cos'l'" - cos6sin$"cos* 
sin 9 sin $" 



(61) 
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cos * sin $" + cos Q cos $" sin sin sin 9 \ 
- sin * sin $" + cos 6 cos $" cos ^I^ cos * sin 8 [62) 
- sin 6 cos $" cos 9 / 



From equations (|5l|), ( ^9| ) and (61), the matrix C is given 
by the product 



C ^A-^B 



(63) 



Performing the right-hand-side product and identifying 
each matrix element by those of expression (62) leads to 
the following trigonometrical relations 

cos 9 = — cos S cos I cos H{t) — sin 5 sin I (64) 

sin 9 cos a = cos 5 sin / cos H{t) — sin 6 cos I (65) 

sin 9 sina = cos Jsini? (t) (66) 

sin 9 sin = — cos ip cos I sin H{t) 

— sin -0 sin S cos Z cos H{t) + sin ■0 cos S sin Z (67) 

sin 9 cos 5* = cos i/j cos 5 sin / 

+ sin -0 cos I sin H{t) — cos ip sin 5 cos / cos H{t), (68) 

where we have used the relation ( ^ ) between $" and a 
and have introduced the local hour angle of the source: 

H{t) = T{t) ^a = at + rGrcc„(0) -L-a . (69) 

From the above relations the beam-pattern factors 
i^+(9,$,*) and Fx (9,$, 5') can be computed for any 
instant t, given the position (/, L) of the detector on the 
Earth, its orientation A, and the position {a,d) of the 
source on the sky, as well as the angle that the vec- 
tor forms with the intersection line of the wave plane 
and the Earth's equatorial plane. First one computes the 
local hour angle H{t) by Eq. (|69|), and the angle 9 by 
Eq. (|64|). The angle $ is computed by means of Eqs. (p5|), 
( |65| ) and (66). Finally, the local polarization angle 'J is 
deduced from Eqs. (67) and (68). 

5. 2. Signal from a rotating magnetized neutron star 

According to Eqs. (38)- (39), the gravitational wave sig- 
nal from a rotating neutron star slightly deformed by its 
magnetic field is 



h+{t) - 


cos fl{t — to) - 




COS 2Q,{t ~ 


to) 


(70) 


hx{t) = 


hxi smCl{t — to) - 


- hx2 


sin2f7(i- 


to) , 


(71) 


with 












h+i = 


~ sin j cost 

ho 

2 tan a 








(72) 


h+2 = 


-ho{l +cos^i)/2 








(73) 


hxi = 


-- sin i 

ha 

2 tan a 








(74) 


hx2 = 


—ho cos i 








(75) 


ho := 


6/3 — p 
cr F 








(76) 



Note that in Eqs. (70)-(71) to accounts for a different time 
origin between the neutron star frame (where Eqs. (38)- 
(39) have been derived) and the detector frame. Formally 
to accounts also for the propagation delay r/c. 

According to the above formulae, when the location 
(a, 6) of the source is known, the signal h{t) measured by 
the detector depends on four a priori unknown parameters: 

— the inclination angle i of the line of sight with respect 
to the rotation axis of the star; 

— the angle a between the magnetic axis and the rotation 
axis; 

— the polarization angle (cf. Figs. and ||) : ip measures 
the orientation of the major axis of the ellipse formed 
by the projection of the neutron star's equator on the 
plane of the sky; 

— the time to such that to — r/c is the instant when 
the magnetic dipole moment vector is in the plane 
formed by the rotation axis and the line of sight, and 
when the scalar product Ad ■ ej is positive. 

The signal measured by the detector can be written with 
the explicit dependence upon these parameters: 

h{t) = [ ,0 ) a) cos f2(t — io) 

+Fx hxi{i,a)smQ{t - to) 



-F+[—,i;] h+2{i)cos2nit~to) 



+Fx yj^'^J hx2{i)sin2n{t-to) (77) 

where tsu is the duration of one sidereal day: and Fx 
are periodic functions of t/tsu^ with period one. 

Let us introduce the amplitudes of the signal at the 
frequencies and 2Q respectively: 



Ai := 



A. 



Fl h\. 



Fx h\i 



Flh\2 



Fx h^x2 



(78) 
(79) 



The daily variation of Ai and A2, computed from 
Eqs. (48)- (49) and (||)-(||), is represented in Figs. |l^ 
for the specific case of the VIRGO detector [L — —10° 30', 
I = -1-43° 40' and A = -161° 33' (P. Hello, private com- 
munication)] and the Crab pulsar (a = 5h34min, S = 
+22° 01'). Note that generally the amplitude of the signal 
at both Q and 2Q is maximum within a few hours of the 
instant when the Crab pulsar crosses the local meridian 
(local sidereal time T = 5 h 34 min) . Figs. p^[]l^ show that 
the measure of the time-varying amplitude at one of the 
two frequencies and 2U allows to determine the polar- 
ization angle i/j and the inclination angle i. If the signal 
is recorded at both frequencies, the angle a between the 
magnetic axis and the rotation axis can be determined as 
well. Note that this angle is a fundamental parameter for 
the theory of pulsar magnetospheres. 
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Fig. 10. Relative amplitude of the signal at the frequency Q {Ai, top) and 2^1 (A-z, bottom) as a function of the local UT time 
t (assuming a vanishing local sidereal time at h UT), for the specific case of the VIRGO detector and the Crab pulsar. The 
polarization angle is V' = 0°. Each plot corresponds to a given value of the inclination angle i. Various lines in the same plot 
correspond to various values for the angle a between the magnetic axis and the rotation axis: solid line: a = 15.00°, dashed 
line: a — 33.75°, dot-dash-dot-dash line: a = 52.50°, dotted line: a = 71.25°, dash-dot-dot-dot line: a = 90.00°. ho is defined 
by Eq. (76). 



6. Conclusion 

We have considered the gravitational radiation emitted by 
a distorted rotating fluid star. The distortion is supposed 
to be symmetric with respect to some axis which docs not 
coincide with the rotation axis. The gravitational emis- 
sion takes place at two frequencies: Q and 20, where Cl 
is the rotation frequency, except in the particular case 
where the distortion axis is perpendicular to the rotation 
axis (only the frequency 20 is then present). As an appli- 
cation, the magnetic field induced deformation is treated. 
If, as usually admitted, the period derivative, P, of pul- 
sars is a measure of their magnetic dipole moment, the 
gravitational wave amplitude can be related to the observ- 
able parameters P and P of the pulsars and to a factor 
p which measures the distortion response of the star to a 
given magnetic dipole moment. (3 depends on the nuclear 



matter equation of state and on the magnetic field distri- 
bution. The amplitude at the frequency 20, expressed in 
terms of P. P and /3, is independent of the angle a be- 
tween the magnetic axis and the rotation axis, whereas at 
the frequency O, the amplitude increases as a decreases. 

Using a numerical code generating self-consistent mod- 
els of magnetized neutron stars within general relativity, 
we have computed the deformation for explicit models of 
the magnetic field distribution and a realistic equation 
of state. It appeared that the distortion at fixed mag- 
netic dipole moment depends very sensitively on the mag- 
netic configuration. The case of a perfect conductor in- 
terior with toroidal electric currents is the less favorable 
one, even if the currents are concentrated in the crust. 
Stochastic magnetic fields (that we modeled by consider- 
ing counter-rotating currents) enhance the deformation by 



14 



S. Bonazzola & E. Gourgoulhon: Gravitational waves from pulsars 



15 



f = 30.00 °, i = 0.00 



f = 30.00 



45.00 



Tp = 30.00 °, i = 90.00 



I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' 







i/ = 30.00 



10 15 20 
t [h] 



0.00 




5 10 15 20 
t [h] 




5 10 15 20 
t [h] 

f = 30.00 °, i = 45.00 




5 10 15 20 
t [h] 




5 10 15 20 
t [h] 

Tp = 30.00 °, i = 90.00 




5 10 15 20 
t [h] 



Fig. 11. Same as Fig. ^ but for the polarization angle = 30°. 



several orders of magnitude and may lead to a detectable 
amplitude for a pulsar like the Crab. As concerns super- 
conducting interiors — the most realistic configuration for 
neutron stars — we have studied numerically type I su- 
perconductors, with a simple magnetic structure outside 
the superconducting region. The distortion factor is then 
~ 10^ to 10'^ higher than in the normal (perfect conduc- 
tor) case, but still insufficient to lead to a positive detec- 
tion by the first generation of kilometric interferometric 
detectors. We have not studied in details the type II su- 
perconductor but have put forward some argument which 
makes it a promising candidate for gravitational wave de- 
tection. Due to the complicated microphysics involved in 
type II superconductors we delay their study to a future 
paper. We also plan to study the deformation induced by 
a possible ferromagnetic solid interior of neutron stars, as 
well as the effects of a strong toroidal internal magnetic 
field. 

Regarding the reception of gravitational waves from a 
pulsar by an interferometric detector, we have computed 
the amplitude modulation of the signal induced by the di- 



urnal rotation of the Earth. By inspecting the wave form, 
and assuming the position of the pulsar to be known (the 
pulsar can be recognized by its period), one can determine 
the inclination angle of the line of sight with respect to the 
pulsar rotation axis, as well as the orientation of the pul- 
sar equatorial plane. Moreoever, by comparing the wave 
forms at the frequencies fl and 2f2, the angle a between 
the rotation axis and the magnetic axis can be determined. 

Pulsars may be good candidates for the detection by 
the forthcoming VIRGO and LIGO interferometric detec- 
tors. A frequently invoked mechanism for gravitational 
emission concerns asymmetries of the neutron star solid 
crust and the resulting precession. In this article, we have 
examined instead the bulk deformation of the star induced 
by its own magnetic field. For some configurations of the 
magnetic field (stochastic distribution, type II supercon- 
ductor), the deformation may be large enough to lead to 
a detectable signal by VIRGO, with the total magnetic 
dipole moment (or equivalently the surface magnetic field) 
keeping its (relatively small) observed value. The positive 
detection of gravitational waves from pulsars would lead 
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Fig. 12. Same as Fig. but for the polarization angle ip = 60°. The case t/j = 90° is identical to ?/> = 0° (Fig. ^ 



to some constraints on the internal magnetic field distri- 
bution, which would be of great interest for the theories of 
pulsar magnetospheres. This would constitute an example 
of a significant contribution of gravitational astronomy to 
classical astrophysics. 
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A. From QI to ACMC coordinates 

Most studies of stationary rotating neutron stars make 
use of quasi-isotropic (QI) coordinates {t,r',9',ip) (cf. the 
discussion in Sect. 2 of Bonazzola et al. 1993). In these 
coordinates, the spatial components of the metric tensor 



have the following asymptotic behavior: 



^ ^ ai ^ Pii cos 29' + /3io 



-3 (Al) 



ge'B' 
9 VP 



m Pii cos26>' + PiQ 



r'' sin^ 9' 



r'^ (A2) 
(A3) 



where ai, as, /3io, and P30 are some constants. By 
comparison with the definition (4), it appears that such 
coordinates are not ACMC to order 1: in order to be so, 
the l/r'^ term in g^v and gg/g' should not contain any 
cos 29', i.e. Pu should vanish. It can be seen easily that the 
following coordinate transformation leads to an ACMC 
coordinate system (<, r, 9, ip): 



Pii cos^ 9 
r 

Pii cos 9 sinf 



(A4) 
(A5) 
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By computing the goo component in the (t, r, 6, ip) coor- 
dinates from the components ga'/S' in the {t,r' ,9' , if) co- 
ordinates and by identification with Eq. (2), one obtains 
the following value of Thome's mass quadrupole moment: 



T 



T- — 



-^yy ~ 2 



A//3ii) 
if z^j, 



(A6) 

(A7) 
(A8) 



where (i) 72 is the coefficient of cos 26*' in the term of 
the 1/r' expansion of the metric component goo in the QI 
coordinates {t,r',9',(p) and (ii) M is half the coefficient 
of 1/r' in the same expansion (M is nothing else than the 
total gravitational mass of the star) . 

To summarize, Thome's quadrupole moment compo- 
nent Jzz can be computed via equation (A7) by reading 
off the coefficients Pn and 72 in the 1/r' expansions of the 
metric components in the QI coordinates. 
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